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Abstract 

A fully non-perturbative lattice determination of the up/down and 
strange quark masses is given for quenched QCD using both, 0{a) im- 
proved Wilson fermions and ordinary Wilson fermions. For the strange 
quark mass with 0(a) improved fermions we obtain m*^'^(/i = 2GeV) = 
105(4) MeV, using the interquark force scale tq. Due to quenching prob- 
lems fits are only possible for quark masses larger than the strange quark 
mass. If we extrapolate our fits to the up/down quark mass we find for 
the average mass mf^{^ = 2 GeV) = 4.4(2) MeV. 
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1 Introduction 



Some of the least known parameters in the Standard Model are the light quark 
masses ruu, md and mg. Their phenomenological values have been discussed since 
the early days of the quark model. Paradoxically, the values of the later discovered 
heavier quarks are more accurately known |l| . The reason is that the connection 
between light quark masses and observables is highly non-perturbative. This 
means that the lattice approach is an appropriate technique for this problem. 

In this paper we shall present a completely non-perturbative determination of 
light quark masses. The recent major step forward has been the non-perturbative 
lattice determination of the renormalisation constants of the mass operators. 
Also, due to the increase in available computer time, a more reliable continuum 
extrapolation is now possible. 

This paper is organised as follows. In section |^ we discuss the definition of 
the quark mass and its renormalisation group behaviour. Transcribing lattice 
data to physical units requires a scale to be set. For quenched QCD this problem 
is discussed in section |[ The lattice technique for obtaining the quark masses 
and their renormalisation is presented in section ^. In section ^ we give our re- 
sults, and in section ^ we extrapolate them to the continuum limit to remove 
residual discretisation effects. We perform the calculations for both, 0{a) im- 
proved fermions and for Wilson fermions. Both should extrapolate to the same 
continuum result, and thus we have a consistency check between the two meth- 
ods. We have previously used tadpole improved perturbation theory to compute 
the renormalisation constants. In section ^ we test the validity of this approach. 
Finally, in section ^ we give our conclusions. 

2 Defining the quark mass 

Due to confinement quarks are not eigenstates of the QCD Hamiltonian and are 
thus not directly observable. A definition of the quark mass from an experiment 
thus means prescribing the measurement procedure. Theoretically this is equiv- 
alent to giving a renormalisation scheme S and scale M. Conventionally, quark 
masses are given in a mass independent scheme, such as the MS scheme, at 
some given scale /x, commonly taken as 2 GeV . In a general mass independent 
scheme 5 the renormalised quark mass is given by 

m^{M) = Zi{M)m,are. (1) 

The running of this renormalised quark mass as the scale M is changed is con- 
trolled by the (3 and 7 functions in the renormalisation group equation. These are 
defined as scale derivatives of the renormalised coupling and mass renormalisation 
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where the bare parameters are held constant. These functions are given pertur- 
batively as power series expansions in the couphng constant. The expansion is 
now known to four loops in the MS scheme B H]. We have 
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with (3 = 1.20206 . . . and (5 = 1.03693 . . ., C being the Riemann zeta function. 
We may immediately integrate eq. (0) to obtain 
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The renormalisation group invariant (RGT) quark mass is defined from the renor- 
malised quark mass as 



m 



RGI 
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where 



[AZg{M)]-' = [26o/(M) 



^mO ( rg^{M) 

" exp 



7m (0 , ^ 
bo^ 



(8) 



(9) 



and the integration constant upon integrating eq. (|^) is given by A*^, and similarly 
from eq. we have m^^^ . A.^ and m^^^ are independent of the scale. Under a 
change of variable (scheme change or iS — iS'), 



G(/) = /(I + ci(/)2 + ...). 



(10) 
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Figure 1: One-, two-, three- and four-loop results for [AZ^^(/i)] ^ in units of 

It can be shown that the first two coefficients of the /3 function, the ffist coefficient 
of the 7 function and m^^^ are independent of the scheme, while A only changes 
as A"^' = A'^exp(ci/&o)- 

For the MS scheme computing [AZ^^(/x)]~^ involves ffist solving eq. (0) for 
g'"^{fi) and then evaluating eq. (|^). This gives Fig. |1|. We expand the j3 and 7 
functions to the appropriate order and then numerically evaluate the integrals. 
At yU = 2 GeV we have fi/A^''^ ~ 8, and it seems that already at this value we 
have a fast converging series in loop orders. Indeed, only going from one loop to 
two loops gives a significant change in [AZ^^(8A*^'^)]~-'^ of order 7%. From two 
loops to three loops we have about 2%. The difference between the three-loop 
and four- loop results is 0(0.5%). So if we are given m^^^ , and we wish to find the 
quark mass in the MS scheme at a certain scale, we need only use the four-loop 
result from eq. @ or equivalently Fig. |l|. 

3 Digression: which scale to use? 

We always need one (or more) experimental numbers as input to set the scale. 
Ideally, it should not matter what quantity we use. Obvious choices are the force 
scale vq [0, or the string tension or some particle mass (e.g. the proton, 
or for quenched QCD at least, the p). So a ffist requirement is that whatever 
quantity we use, we should be in a region where the scaling to the continuum 



4 



1.25 



1.20 - 



1.15 - 



1.10 - 



1.05 ^ 
0.85 



(rya) X (aVo) 



0.90 



0.95 



1.00 



1.05 



1.10 



Figure 2: The product of r^/a and the string tension ay/a. {ro/a){gQ) is taken 
from the formula given in J^, while the string tension is taken from ((3 = 5.7 , 
5.8, 5.9, 6.4;, ^ ([3= 6.0, 6.2;, gj (13 = Q.A), ^ ((3 = Q.Q). 



limit is the same for all quantities. Thus for tq and the string tension we wish 
that 

— {go) X {a^/a){go) = const. (11) 
a 

over the g^ = 6//3 region used in the simulations, and indeed for all smaller g^. 
(We know that this must break down below a value of P around 5.7, due to the 
appearance of non-universal terms.) In Fig. § we show this product. This seems 
reasonably constant, with a fit value of 1.170(5). 

The second requirement is to set the scale in MeV. As we are considering 
quenched QCD, it is not obvious that choosing scales from different experimental 
(or phenomenological) quantities will necessarily lead to the same results. Indeed, 
in the real world 0, the values are 

ro = 0.5fm= (394.6 MeV)-\ 

= 427 MeV, (12) 

(lfm~^ = 197.3 MeV) which gives for the product a value of 1.082 - almost a 10% 
difference from the quenched value. As both phenomenological estimates come 
from the same potential model [lT]|, presumably the quenched lattice potential 



has a slightly different shape from the continuum potential. 
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Recently, the ALPHA collaboration has determined a value for A 



|T2| 



of 



A*^^ = 0.602(48)/ro, (13) 

which may easily be converted using eq. (|TT]) to the string tension scale. However, 
the numerical value will then suffer from the same 10% ambiguity. Thus we find 

= 238(19) MeV, 

= 220(18) MeV. (14) 

Our recent spectrum results |T^, using 0{a) improved fermions, also show a 
difference whether one uses the p or the proton mass to set the scale. Using tq 
to set the scale is roughly equivalent to using m^. 

In the following we shall adopt the tq scale as given in namely 

ln(a/ro) = -1.6805 - 1.7139(/? - 6) + 0.8155(/3 - 6)^ - 0.6667(/5 - 6)^ (15) 

(with an error of 0.3% increasing to 0.6% for f3 in the range 5.7 < (3 < 6.57), 
but delay using a numerical value for this scale for as long as possible. For the 
standard scale of /i = 2 GeV this gives, upon solving eqs. (0) and (H) to the 
appropriate loop order, the results for [AZ^^(/i)]~^ shown in Table [l|. For later 
reference the results for some other fi values are also given there, together with 
af^ip). 





one-loop 


two-loop three-loop 


four-loop 






2.00 GeV 

2.12 GeV (1/a at p = 6.0) 
2.90 GeV (1/a at (3 = 6.2) 
3.85 GeV (1/a at p = 6.4) 


0.760(10) 
0.752(10) 
0.716(8) 
0.689(7) 


0.704(9) 
0.697(8) 
0.667(7) 
0.644(6) 


0.718(10) 
0.711(9) 
0.677(7) 
0.652(5) 


0.721(10) 
0.714(10) 
0.679(8) 
0.653(6) 




af-(/i) 


2.00 GeV 

2.12 GeV (1/a at p = 6.0) 
2.90 GeV (1/a at p = 6.2) 
3.85 GeV (1/a at P = 6.4) 


0.268(10) 
0.261(9) 
0.228(7) 
0.205(6) 


0.195(6) 
0.191(5) 
0.170(5) 
0.156(3) 


0.201(6) 
0.196(6) 
0.174(5) 
0.159(4) 


0.202(7) 
0.197(6) 
0.175(5) 
0.159(4) 



A MS 

A^ 



Table 1: Useful values of [AZ^^'ifi)]-^ and af^fi) = {g"%fi)y/An. The errors 
are a reflection of the error in eq. (|7^. The values ofl/a are found from eq. ( |7^ 
together with tq from eq. (|7^j. 
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4 Determining the quark mass 




We shall now derive formulae for the quark masses using the conserved vector 
current (CVC) and the partially conserved conserved axial vector current (PCAC) 
by assuming Taylor expansions in the bare quark mass for the relevant functions 
that occur. We distinguish two quark masses. The Ward identities arising from 
an infinitesimal vector transformation in the partition function lead to a bare 
quark mass given by 

^ = 1,2, (16) 

where Kg^ is the corresponding hopping parameter, and is the critical hopping 
parameter. This is the standard definition of the quark mass. Similarly, for an 
infinitesimal axial transformation the Ward Identity (WT) or PCAC definition of 
the quark mass can be written as 

~ , ~ t»0 (a4^f'^^(t)P''l^^'— ^'^(0)) ,^ , 

amo^ + arria^ = —— -— — ,, ,. (17) 

= — smh arupq , 18 

{0\V'^i<i2\pS) ^ ' 

where A (V) is the axial vector current (pseudoscalar density). The precise form 
of A and and V will be given later for the 0{a) improved as well as the Wilson 
cases. We have summed the operators over their spatial planes. While A{t) and 
V(t) should be point operators, to improve the signal P(0) is smeared over its 
spatial plane. To obtain the second equation, we have re-written eq. (^) in a 
Fock space and then introduced a complete set of states in the usual way. We 
have then picked out the lowest pseudoscalar (PS') or 0^+ state whose mass we 
denote by rripg^. 

Both definitions of the quark mass must be renormalised. In a scheme S at 
scale M we have 

ml{M)+ml{M) = Zi{M,amg,,amg,) {rfig, +mg,) . (19) 

The Ward identities give = '^/Zs (from CVC) and Z^ = Za/Zp (from 
PCAC). 

Let us now Taylor expand m and the pseudoscalar mass mp|^ in terms of the 
bare quark masses m^-. 



Y 



1 + am„, + arrio, + a- — , — + . . . 



\{amq^ + arriq^), 



7 



[am 



<?1'?2N2 

PS ) 



Y, 



PS 
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Maniq^ + am-gj. (20) 



The functions must be symmetric under interchange of the quarks, i.e. qi ^ q2- 
Only at the lowest (first) order in the quark mass is the functional form simply 
anig-^ +amq^. At the next order both terms, (am^^ +amq^Y {aniq^Y + {aniq^)'^ 
are allowed. Taylor expanding eq. (|19D and comparing with eq. (pOD gives us 



Y 



Zi{M) ZI{M)Za 



(21) 



Renormalisation constants which show no explicit quark mass dependence refer 
to rUq^ = rUq^ = 0. 

We shall also Taylor-expand the matrix elements appearing in eq. (0). First 

and 
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we define the bare pseudoscalar decay constant by {0\Af^\PS) = m 
similarly we set {0\V'^^'^'^\PS) = —gps^- Expanding the decay constants fps'^ and 
Qpg^ to first order in the quark masses gives 
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(22) 



(23) 



Thus, at least to this order, we have a relation between the WI quark masses and 
the pseudoscalar mass. 

Using eq. ([19|) gives the renormalised quark mass, and we additionally use 
eq. (13) to re-write this in a RGI form as 
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Upon taking the continuum limit Qq 0, any scaling violations will show them- 
selves as non-constant Ca, Cb functions. Equation ( P^ is the main result of this 
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analysis. Given c*, cl and the pseudoscalar mass, we can then determine the 
quark masses. 

For the {us) we set qi = u, q2 = s, and for the K'^ (ds) we set qi = d, 
q2 = s. Together with the (ud), with qi = u, q2 = d, this gives from eq. (|2^) 



rom 



RGI 
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+ ... 



RGI 



(27) 



where we have defined mf'^^ 



[m. 



RGI 



_l_ uiR<^^y2, i.e. the average of the u/d 



quarks. We have ignored any small corrections due to electromagnetic effects. 



5 Numerical results 



5.1 Pseudoscalar mesons and bare quark masses 



For degenerate quark masses from eq. ([20|) we have 



aruq = Y [1 + {c + d)amq + . . .] am,, 
{anipsY = Yps [1 + {cps + d)amq + . . .] am, 



(2^ 



and 



where mps 



anin 
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{ampsY Yps 



l+{^^]{ampsr + ... 



Y, 



PS 



(29) 



m 



QQ 



ipg [I.e. qi = q2 = q). Equation (|^) gives Y /Yps for the Ca term 
in eq. (P?]), but the gradient (c — Cps)/Yps is not sufficient to give —Cps/Yps 
for the Cfe term. For 0(a) improved fermions, associating the mass expansion 
parameters Ba, bp and bm with our expansions, we find c = — (6^ — bp) and 
d = bm. First order perturbation theory |]T5| gives c ~ 0.001(y'Q. Q On top of that 
c/cps = 0{a), so that the effect of c can safely be ignored. For Wilson fermions 
we shall assume that either c is small in comparison with cps, as above, or that 
the complete term Cblromx)'^ is small when compared with c^. As we shall see, 
little error is introduced by this assumption. 

We have computed the pseudoscalar mass nips and the Wbare quark mass 
both, for 0(a) improved fermions and Wilson fermions. For improved fermions 
the calculations were done at /3 = 6.0, 6.2, 6.4 and csw = 1.769, 1.614, 1.526 [0, 
respectively, while for Wilson fermions we only did calculations at /? = 6.0, 6.2. 
The computational methods used are standard. For the pseudoscalar mass we 



non-perturbative estimate |16 , however, gives c ~ —0.15 a,t P — 6.2. 
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Figure 3: [ampsY /amq and amg/arriq against am q for O [a) improved fermions. 
Filled circles denote points used in the fits. 



used the correlation function 



i>0 



A 



-nipst _|_ ^-mps{T-t) 



(30) 



T being the temporal extent of the lattice. To improve the signal, a Jacobi- 
smeared operator was used, as described in [|T0|. For Wilson fermions the pseu- 
doscalar meson mass results can also be found in [|l^]. The 0(a) improved results 
for mps are given in Table ^ in the Appendix. The various k values used, the 
lattice size, and the number of configurations generated are also collated there. 

For aniq the ratio of two point correlation functions as given in eq. ([T^ was 
used. For 0(a) improved fermions, as well as the action, the operators must also 
be improved: 

Af, = + CAad^P, 
V = P, (31) 

where = ^7^75*3' and P = 5755. By choosing the improvement coefficient 
CAido) appropriately, the Ward identity can be made exact to 0(a). CA{go) is 



non-perturbatively known |jT^. In Table y in the Appendix we give our results 
for aniq. 

Let us first discuss 0(a) improved fermions. In Fig. ^ we show the ratios 
{ampsY / aniq and amq/arriq against arriq, while in Fig. ^ we plot arfiq/ {ampsY 
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Figure 4: amq/ {ampsY against ((ro/a) x {amps)Y for 0{a) improved fermions. 
Filled circles denote points used in the fits. The dashed line 3.13j is mps = 
\/2mK (which here corresponds to a fictitious ss bound state) while the dotted 
line 0.125j is mT^. 



against ((ro/a) x [amps)Y. We must now search for a region where eq. 
without higher order terms, is vahd. For large quark mass values we expect 
non-linear terms, while for small quark masses quenched QCD chiral logarithms 
become significant. Finite volume effects do not seem to be a problem, as for 
(3 = 6.0, K = 0.1342 and (3 = 6.2, k, = 0.1352 we have made runs on two different 
volumes, without significant changes in the results. 

We now make some cuts. In Fig. ^ we see that for small quark masses there 
are significant deviations from linearity. In particular the light quark mass mi 
lies in a region where no direct linear extrapolation is possible. However, above 
nips ^ \/2mK deviations from linearity seem small. We shall thus assume that 
at least above the strange quark mass any effects of chiral logarithms are small. 
For heavy quark masses, on the other hand, linearity is still present until at least 
rUq ~ Snis ~ ^rric. (Note that 2(rom£))^ ~ 44.9). In this interval lie four or more 
quark masses. The results of the various fits are given in Table ^. 

For Wilson fermions we simply set = A^. The results for aniq have also 
been given in |T0|. In Figs. |] and || we show the ratios (arrips)^ / aniq, arhq/amq 
and amq/ {anipsY ■ Similar fit ranges as for improved fermions seem appropriate 
once more. 

To illustrate the dependence of some of these results, we show in Fig. ^ 
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Y 
Yps 


Y ^. c-cps (ro\ ^ 
Yps Yps \aj 


Y 


0(a) improved fermions 


6.0 
6.2 
6.4 


0.314(2) 
0.465(2) 
0.617(5) 


-0.0017(2) 
-0.0042(2) 
-0.0063(5) 


0.972(6) 
1.060(3) 
1.082(5) 


Wilson fermions 


6.0 
6.2 


0.368(4) 
0.502(6) 


-0.00025(51) 
-0.0036(5) 


0.711(12) 
0.814(10) 



Table 2: Fit values. 
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Figure 6: The same as Fig. El hut for Wilson fermions. 



the results for Y taken from Table | 
form 



I together with Pade-like interpolations of the 

1 + Pi9o + P29o 



1 + (pi - c)g^ + psg^ 



4' 



(32) 



arranged so that the perturbative result Y{go) = 1 + cgQ with c = 0.09051 for 
improved fermions and 0.05195 for Wilson fermions is obtained for small g^. 
Possible Fade interpolations are found to be (pi,P2,P3) = (—1.24,0.256,0.347) 
for 0{a) improved fermions and (—0.944,0.00,0.0746) for Wilson fermions. Also 
shown for comparison are 0{a) improved results found in |16|. We see that for 



0{a) improved fermions first order perturbation theory is good for gQ < 0.96, 
while for Wilson fermions a breakdown occurs much earlier. 

Thus we now have estimates for Y/Yps and (c — Cps)/Yps- Y will also be 
needed for Wilson fermions. 



5.2 Renormalisation 

To compute and we must now determine the factor AZ^{M)Z^{M). For 
0{a) improved fermions this was done by the ALFHA collaboration [|12| using the 
Schrodinger Functional (SF) method. With the notation of eq. (§) their result 
can be written as 

AZ^^{L~^)Z^^{L~^) = 1.752 + 0.321(/? - 6) - 0.220(/? - 6)^ (33) 
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Figure 7: Y against Qq. Our 0{a) improved fermion results are shown as filled 
circles, while those from pm are shown as open circles. The Wilson fermion 
results are filled squares. The one-loop perturbation theory results are also shown 
for the 0{a) improved case (dashed line) and the Wilson case (long dashed line). 



(valid for 6.0 < (3 < 6.5), where they have worked at a scale given in terms of 
the box size L. As emphasised previously, this is a mapping from the bare quark 
mass to the RGI mass, so this function is the same in all schemes. In eq. (0) 
the total error is about 2%. 



For Wilson fermions we use the method proposed in [|1^] and refined in 
This mimics perturbation theory in a certain MOM scheme by considering am- 
putated quark Green's functions in, say, the Landau gauge, with an appropriate 
operator insertion. The renormalisation constant is fixed at some scale p^. This 
gives a non-perturbative determination of Zg°^'{p). {Zp°^'{p) is not suitable, 
as chiral symmetry breaking means that Zp°^''- ^ as we approach the chiral 
limit, as recently emphasised in pOf-) More details of the method, our momen- 



tum source approach, and results are given in [|T9|. As AZ^^ is known in the 
MS scheme at scale /i (see Fig. ^ and Table |I]), we can write 

^zT'\v)zi°'\v) = -^=r§ ^ ^(/^)' (34) 

where we have used eq. (^) and the definitions = '^/Zs, AZ^ = 1/AZs and 

zri^^)=xs.;^oM{f^,p)zrip), (35) 
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with 



^S:MOMif'^P) 



(36) 



Here X converts the renormahsation constant from the MOM scheme to the 
MS scheme and can be calculated using a continuum regularisation (e.g. naive 
dimensional regularisation). So we can write 



7 MOM 



1 + 



a" 



47r 



a' 



An 



B2 + ... 



(37) 



with Bi = 16/3 |]To|- The coefficient B2 has recently been calculated [21|, giving 
in the Landau gauge for Nf = flavours a value of 177.48452. Hence, knowing 
the coefficients Bi, B2, we can trace back to the more general expression for X. 
A suitable formula is 



[^S:MOMif^^P)] ^ 



exp 



^(0 



MS 

(]f Jul— 111. -L 



7^(0 - j(°"(g(0) 



(3^ 



^gMOM _ Q(^gMs^ from eq. (|10D)- Armed with this estimate for X, then from 
eq. ( P^D we see that Zg^{fi) should be independent of p^, so plotting Zg^{n) 
against we expect to see a plateau. In Fig. |] we show this, plotting first the 
original data, extrapolated to the chiral limit, then the results for Zg''^ when only 
using Bi, and finally using both Bi and B2. We see the data for Zs becoming 
flatter. The results of the fit are given in the ffist column of Table |^. The 
appropriate values for Y (see Table 0) and AZ^^(/i = 1/a) (see Table ^ are 
substituted in eq. ( p4|) to give the results listed in the second column of this 
table. 



p 


Zr(/i = l/a) 


AZ*--(p)Z-°-(p) 


6.0 

6.2 


0.790(5) 
0.803(5) 


2.49(6) 
2.25(4) 



Table 3: Fit values for the Wilson renormalisation constants. 



6 Continuum results 



Plotting Ca and Ch against a? for 0(a) improved fermions, and against a for Wilson 
fermions, we can now extrapolate these coefficients to the continuum limit. In 
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0.9 
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zr{^=i/a)i 



1 loop 



- , -5- -n ' 



-.MOM, , 

Zs (P) 
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10.0 
p' [GeV'] 



100.0 



Figure 8: Zg^{ii = 1/a) against for /5 = 6.0. (See Pm for details.) The open 
squares are the original data, in the chiral limit, while the filled circles represent 
the results of multiplying by X using both coefficients Bi and B2- The 

open circles are the result of using only the Bi coefficient. The straight line is a 
fit to the plateau, the length denoting the fit range chosen. 
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Figure 9: The continuum extrapolation for Ca and Cb for 0{a) improved fermions. 
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Figure 10: The continuum extrapolation for Ca and Ch for Wilson fermions. 



Fig. H we show the results for Ca and Ch for improved fermions. A hnear fit is also 
plotted. The results of this fit are given in Table ^. As anticipated, using the 
first order perturbative result from for c in Cf, has no influence on the result. 
In Fig. [TO we show the equivalent results for Wilson fermions. In this case, as 
we only have two (3 values, the fit degenerates to an extrapolation. The results 
of this extrapolation are also given in Table |^. We note that the results for 0(a) 
improved fermions and Wilson fermions for c* are compatible with each other. 





c* 


r* 


0(a) improved 
Wilson 


0.124(2) 
0.107(17) 


-0.0015(1) 
-0.0038(9) 



Table 4: The continuum extrapolation of Ca c* and Cb cl for 0(a) improved 
fermions and Wilson fermions. 



Upon inserting these numbers in eq. (|27|), we find our estimate for the RGI 
strange and u/d quark masses. To convert to physical numbers, we now have 
in quenched QCD the uncertainty in the scale, as discussed in section 0. If we 



use the scale tq = 0.5 fm, eq. (p!2D, then together with the experimental values 
of the 71 and K masses, namely mT,+ = 139.6 MeV and mK+ = 493.7 MeV, 
rrixo = 497.7 MeV, we find the results for 0(a) improved fermions 



RGI 

s 

mf^^ = 6.1(2) MeV. (39) 



^HGi ^ 146(4) MeV 

,RGI 



^Note, however, that the non-perturbative estimate from ||T^ at /3 = 6.2 is c/Yps x 
(ro/a)~^ ^ —0.0012 and has a similar order of magnitude to the kept term. 
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The error comes from Table ^ and from eq. (|33| ) . As can also be seen from Table ^ 
most of the result comes from the constant term, with the slope giving only a 
small correction to the answer. 

For Wilson fermions we have mf"^^ = 121(20) MeV, mf^^ = 5.3(8) MeV. 
This result is somewhat lower than the 0{a) improved numbers. We ascribe this 
mainly to the fact that we only have two values of (3, which makes a continuum 
extrapolation more difficult. Also the number of n values used and the size of 
the data sets are smaller than for 0{a) improved fermions. Nevertheless, within 
a one-standard deviation the results are in agreement. 

In the MS scheme at the 'standard' value of /x = 2 GeV, using the four-loop 
results from Table |^, we find for 0(a) improved fermions 

mf^(/i = 2GeV) = 105(4) MeV, 

mp(/i = 2 GeV) = 4.4(2) MeV. (40) 

The corresponding Wilson results are mf^{fi = 2 GeV) = 87(15) MeV and 
mp(/i = 2 GeV) = 3.8(6) MeV. 

Note that for the mi quark mass result we have simply extrapolated the fits 
for the strange quark mass result downwards. The mass ratio ms/mi for 0{a) 
improved fermions is ~ 23.9, which is very close to the value given in leading 
order chiral perturbation theory - namely {rn\-+ + mj^o — 'm'^+)/m'^+ ^ 24.2 (see 
eq. (l27|)). This is simply because |c^| <^ |c*|/(romx)^, and so the second term in 



eq. ([27|) is almost negligible. The mass ratio is then independent of c*. 



In Fig. 11 we plot our results. Below the dotted line, we have given our 



previous result using the string tension as the scale, as given in eq. ([l^ . 
As a comparison we have also re-plotted our result given in eq. (|iO| ) using ^/a as 
the scale. A reasonable agreement is seen. As our previous result used tadpole 
improved (Ti) perturbation theory to determine the renormalisation constant, 
it would seem that the use of TI perturbation theory does not introduce much 
error. In the next section we shall briefiy investigate this point. 



7 Digression: comparison with tadpole-improv- 
ed perturbation theory 

In this section we shall discuss how reliable TI perturbation theory is. Lowest 
order perturbation theory gives 



Zr(/i = l/a) = l 



IGtt' 



■CFB''%csw) + 0{g',). 



(41) 



The B coefficient is uncomfortably large [|T^. In ||22| this was traced to large 
tadpole diagrams in the perturbation expansion and also to the expansion in 
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Figure 11: The strange and light quark masses from eq. ( Q^ ) ('This work'). The 
Wilson fermion results are shown dotted. Also shown is our previous result ('pre- 
vious') The arrows denote the result when using as a scale. 



a non-physical (bare) coupling constant. Removing the tapole diagrams and 
expanding in (say) af'"^ gives the improved series 



ryMS 



(/i = 1/a) = Mo 1 



a: 



TC 



+ 0{ai), (42) 



with uq 



{^TiUpiaq)'^ to be numerically determined, and csw 



CswU^ 



0- 



A 



description of our variation of this procedure is given in |TD[. In particular, we 
choose Csw to be the non-pert urbatively determined result, rather than setting 
it to be equal to I/mq. af^ corresponds to the four-loop results (see Table |I|) 
also used in . 

As we now have a genuine non-perturbative determination of Zm available, it is 
of interest to compare how the different results scale to the continuum limit. It is 
convenient to first define for two 0{a) improved estimations of the renormalisation 
constant the ratio 

AZi^{M,)Zi^{M,] 



R 



l + 0{a' 



(43) 



AZ^^^{M,)Zi^{M,) 

For ^2 we choose the SF scheme and the result given in eq. (PB|), while for Si we 



use the MS scheme, together with either the perturbative result, eq. (|41|), or the 
TI result, eq. 
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Figure 12: The results for the ratio R as defined in eq. 1^4^) ■ has been 

determined from eq. or eq. (^) and ^Z^^ taken from Table [|. A simple 

interpolation has been used between the f3 values. The dot-dashed vertical lines 
correspond to (3 = 6.0, 6.2, 6.4, from right to left. 



In Fig. |T2| we show the results for the ratio R for perturbation theory and TI 
perturbation theory, using consistently the one-, two- and four-loop results from 
Table [l|. As we originally used one- loop perturbation theory results, it seems more 
consistent to also use the one-loop result for AZ^ when converting Z^^^ to the 
RGI result. This gives the solid line in the figure. This was the approach adopted 
in [|T0|. We expect 0(a) effects to become apparent as — > if Si is not exactly 
0{a) improved. However, a linear fit in for the T/ result (with one-loop AZ^) 
appears to go to i? = 1 with an error of only about 2%, while for the equivalent 
perturbative result the error is about 10%. Thus, in this case tadpole improving 
the perturbative result does give better results. However, choosing other loop 
orders changes the picture somewhat and can make using perturbation theory a 
better choice. We would hke to emphasise that this picture does not have to hold 
for other renormalisation constants. Strictly speaking a case-by-case analysis is 
required. 



8 Conclusions 

In this article we have calculated the strange and u/d quark masses for quenched 
QCD, both for 0(a) improved fermions and Wilson fermions, using a non-per- 
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turbatively determined renormalisation constant. Our results are given in eq. (^) 
and the lines that follow it. Corrections to leading order chiral perturbation 
theory are small, if we stay away from the region where chiral logarithms become 
significant. We have also seen that using TI perturbation theory rather than 
simple perturbation theory does not automatically lead to an improvement of 
the continuum result. 
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Note added 



While this work was being completed, we received a copy of [^. This contains 
some similar results to ours. 
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Appendix 



In Table ^ we give our parameter values used in the 0{a) improved fermion 
simulations together with the measured pseudoscalar mass. For most of the 



overlapping values with [10| there has been some increase in statistics. 
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Volume 


# configs. 


amps 


6.0 


1.769 


0.1217 


16^ 


X 32 


0(150) 


1.2546(12) 


6.0 


1.769 


0.1263 


16^ 


X 32 


0(150) 


0.9704(11) 


6.0 


1.769 


0.1285 


16^ 


X 32 


0(150) 


0.8189(11) 


6.0 


1.769 


0.1300 


163 


X 32 


0(160) 


0.7071(16) 


6.0 


1.769 


0.1310 


163 


X 32 


0(160) 


0.6268(16) 


6.0 


1.769 


0.1324 


163 


X 32 


0(990) 


0.5042(7) 


6.0 


1.769 


0.1333 


163 


X 32 


0(990) 


0.4122(9) 


6.0 


1.769 


0.1338 


163 


X 32 


0(520) 


0.3549(12) 


6.0 


1.769 


0.1342 


163 


X 32 


0(1300) 


0.3012(10) 


6.0 


1.769 


0.1342 


243 


X 32 


0(200) 


0.3017(13) 


6.0 


1.769 


0.1346 


243 


X 32 


0(200) 


0.2390(12) 


6.0 


1.769 


0.1348 


243 


X 32 


0(200) 


0.1978(16) 


6.2 


1.614 


0.1247 


243 


X 48 


O(IOO) 


1.0284(9) 


6.2 


1.614 


0.1294 


243 


X 48 


O(IOO) 


0.7217(9) 


6.2 


1.614 


0.1310 


243 


X 48 


O(IOO) 


0.6043(9) 


6.2 


1.614 


0.1321 


243 


X 48 


0(260) 


0.5183(6) 


6.2 


1.614 


0.1333 


243 


X 48 


0(560) 


0.4136(6) 


6.2 


1.614 


0.1344 


243 


X 48 


0(560) 


0.3034(6) 


6.2 


1.614 


0.1349 


243 


X 48 


0(560) 


0.2431(7) 


6.2 


1.614 


0.1352 


243 


X 48 


0(260) 


0.2016(10) 


6.2 


1.614 


0.1352 


323 


X 64 


O(llO) 


0.2005(9) 


6.2 


1.614 


0.1354 


323 


X 64 


0(290) 


0.1657(6) 


6.2 


1.614 


0.13555 


323 


X 64 


0(280) 


0.1339(7) 


6.4 


1.526 


0.1313 


323 


X 48 


O(IOO) 


0.5305(9) 


6.4 


1.526 


0.1323 


323 


X 48 


O(IOO) 


0.4522(10) 


6.4 


1.526 


0.1330 


323 


X 48 


O(IOO) 


0.3935(12) 


6.4 


1.526 


0.1338 


323 X 48 


0(200) 


0.3213(8) 


6.4 


1.526 


0.1346 


323 


X 48 


0(200) 


0.2402(8) 


6.4 


1.526 


0.1350 


323 


X 48 


0(200) 


0.1923(9) 


6.4 


1.526 


0.1353 


323 


X 64 


0(260) 


0.1507(8) 



Table 5: Parameter values used in the simulations, together with the measured 
pseudoscalar mass. 

The results for the 1^/ quark mass, am^, are first split into two pieces, 
denotes the mass coming from the {diA4^P^^^'^'^^'^) / l^p p^^^^o-^^^d-'^ ratio, while 2am^^) 
is the result of (ylp p^rneared^ ^^p psmearedy r^^^ 2am, = 2amf^ + 2cAamW 
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gives the WI quark mass. All these results are given in Table 0. We define 



/xy 



+4,y 



)/2. 8484 has been replaced by (V; 

2 



2^ 



(^4; 

^x-i,y continuum limit both 8484 and V4 give the same derivative. On the 

lattice we choose the discretisation V| with the smallest (temporal) extension. 
In |T^ the choice 8484 was used. 



/3 


Csw 




2am^^ 

y 


y 


2amg 


6.0 


1.769 


0.1217 


0.9677(7) 


1.7867(33) 


0.8197(5) 


6.0 


1.769 


0.1263 


0.5936(4) 


1.0173(24) 


0.5093(3) 


6.0 


1.769 


0.1285 


0.4340(3) 


0.7080(22) 


0.3754(3) 


6.0 


1.769 


0.1300 


0.3314(3) 


0.5227(21) 


0.2881(3) 


6.0 


1.769 


0.1310 


0.2651(3) 


0.4072(20) 


0.2313(3) 


6.0 


1.769 


0.1324 


0.1751(1) 


0.2609(8) 


0.1535(1) 


6.0 


1.769 


0.1333 


0.1186(1) 


0.1737(7) 


0.1042(1) 


6.0 


1.769 


0.1338 


0.08734(21) 


0.1274(9) 


0.07678(18) 


6.0 


1.769 


0.1342 


0.06282(16) 


0.09217(62) 


0.05519(15) 


6.0 


1.769 


0.1342 


0.06294(27) 


0.0932(12) 


0.05522(25) 


6.0 


1.769 


0.1346 


0.03772(32) 


0.0584(11) 


0.03289(32) 


6.0 


1.769 


0.1348 


0.02485(37) 


0.0399(12) 


0.02154(37) 


6.2 


1.614 


0.1247 


0.7342(3) 


1.1520(18) 


0.6915(2) 


6.2 


1.614 


0.1294 


0.4011(2) 


0.5451(12) 


0.3809(2) 


6.2 


1.614 


0.1310 


0.2966(2) 


0.3776(10) 


0.2826(2) 


6.2 


1.614 


0.1321 


0.2272(1) 


0.2748(7) 


0.2170(1) 


6.2 


1.614 


0.1333 


0.15265(5) 


0.1741(4) 


0.14620(5) 


6.2 


1.614 


0.1344 


0.08542(6) 


0.09336(36) 


0.08196(5) 


6.2 


1.614 


0.1349 


0.05510(6) 


0.05985(34) 


0.05288(6) 


6.2 


1.614 


0.1352 


0.03688(11) 


0.04073(46) 


0.03537(11) 


6.2 


1.614 


0.1352 


0.03703(12) 


0.04046(37) 


0.03553(12) 


6.2 


1.614 


0.1354 


0.02478(6) 


0.02773(18) 


0.02375(6) 


6.2 


1.614 


0.13555 


0.01549(8) 


0.01802(20) 


0.01482(8) 


6.4 


1.526 


0.1313 


0.2709(1) 


0.2870(9) 


0.2637(1) 


6.4 


1.526 


0.1323 


0.2090(1) 


0.2074(9) 


0.2038(1) 


6.4 


1.526 


0.1330 


0.1660(1) 


0.1567(8) 


0.1621(1) 


6.4 


1.526 


0.1338 


0.1172(1) 


0.1045(4) 


0.11461(6) 


6.4 


1.526 


0.1346 


0.06882(6) 


0.05799(38) 


0.06736(6) 


6.4 


1.526 


0.1350 


0.04467(7) 


0.03739(35) 


0.04373(6) 


6.4 


1.526 


0.1353 


0.02661(5) 


0.02276(22) 


0.02603(5) 



Table 6: Results for the Ward Identity quark mass 2amg. 
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